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We study generation of acoustic waves by a processing spin of a nanoparticle deposited on the 
surface of a solid. Our approach elucidates macroscopic dynamics of the Einstein - de Haas effect. 
It is based upon solution of parameter-free equations of motion mandated by the conservation of the 
total angular momentum: spin -|- mechanical. We show that the amplitude of the acoustic waves 
generated by the spin precession may be sufficient for detection by a tunneling electron microscope. 


I. INTRODUCTION 

One hundred years ago Einstein and de Haas had 
demonstrated experimentally that the change in the mag¬ 
netic moment of a freely suspended macroscopic body 
generated its mechanical rotation due to conservation of 
the total angular momentunii. (Back then the concept 
of spin had not existed yet, magnetic moments of atoms 
were attributed to the atomic orbital currents carrying 
angular momentum.) Initially the effect has been inten¬ 
sively studied as a tool to measure the gyromagnetic ratio 
of solids. Subsequent discovery of the electron spin res¬ 
onance (ESR) and the ferromagnetic resonance (FMR) 
provided much more accurate determination of the gyro- 
magnetic ratio, greatly diminishing the significance of the 
mechanical studies based upon the Einstein - de Haas ef¬ 
fect. Meanwhile microscopic mechanisms of the transfer 
of the spin angular momentum to the rotation of the me¬ 
chanical system as a whole remained poorly understood, 
resulting in a number of revisions in recent years^. 

The problem somewhat simplifies when the number of 
mechanical degrees of freedom is finite. This prompted 
a significant body of work on spins in mechanical res¬ 
onators inspired by progress in manufacturing micro- and 
nanoelectromechanical devices (MEMS and NEMS)^. 
Einstein - de Haas effect at the nanoscale has been ex¬ 
perimentally studied in magnetic microcantilevers^ and 
theoretically explained by the motion of domain wallsi. 
Switching of magnetic moments by mechanical torques 
has been proposed^riS. Mechanical resonators containing 
single magnetic molecules have been studied^ " i-—, Exper¬ 
iments have progressed to the measurement of the me¬ 
chanical coupling of a single molecular spin with a carbon 
nanotubei^di. 

The aim of this paper is to generalize the angular- 
momentum-based approach developed for resonators to 
systems with infinite number of degrees of freedom. As 
an example of such problem we study generation of acous¬ 
tic waves by a magnetic nanoparticle on a solid surface 
in the ESR setup. The precession of the spin results 
in the time-dependent mechanical torque acting on the 
nanoparticle. The parameter-free nature of the torque 
can be traced to the conservation of the total angular 
momentum (spin + mechanical). The resulting mechan¬ 
ical vibration of the particle on the surface of the solid 


generates acoustic waves. The effect is proportional to 
the ratio of the magnetic anisotropy energy and elastic 
energy, which has relativistic smallness due to spin-orbit 
origin of the magnetic anisotropy. Nevertheless, the dis¬ 
placement of the surface in the vicinity of the particle is 
sufficient to be detected by the tunneling electron micro¬ 
scope (TEM). 

The following argument makes the acoustic TEM de¬ 
tection of spin precession possible even in the presence 
of stronger effects on the tunneling current. The cur¬ 
rent, being the polar vector that changes sign due to 
spatial reflection, can only depend on the magnetic mo¬ 
ment squared, because the magnetic moment is an axial 
vector that does not change sign under spatial reflection. 
Consequently, the direct interaction between the mag¬ 
netic moment and the tunneling current can only pro¬ 
duce oscillations of the current at double the precession 
frequency. However, generation of the acoustic waves at 
the precession frequency will make the tunneling current 
through the oscillating surface to oscillate at the same 
frequency. Oscillations of the tunneling current at the 
precession frequency in the vicinity of the magnetic clus¬ 
ter has been detected in the pastier—. The proposed 
explanation was based upon oscillation of the tunneling 
conductance due to the spin-orbit coupling between the 
magnetic atom and the tip of the TEM^“— . 

The paper is structured as follows. The model of a 
uniaxial magnetic particle coupled to the elastic matrix 
via the orientation of its easy magnetization axis is de¬ 
scribed in Section [III Bulk waves generated by the vibra¬ 
tion of the nanoparticle due to the precession of its spin 
induced by the ac magnetic field are studied in Section 
EH The boundary conditions and the surface waves are 
worked out in Section lYl Results are discussed in Sec¬ 
tion El where numerical estimates are made for an FePt 
nanoparticle on a gold surface. 


II. MODEL 

The nanoparticle is described within the model of an 
uniaxial magnetic anisotropy, that is, its Hamiltonian op¬ 
erator is given by 

'Hcf = —D{n ■ ^ , 


( 1 ) 
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where n and S represent the anisotropy axis and the vec¬ 
tor of spin operators of the nanoparticle, respectively, 
and D is the anisotropy constant. This spin is embedded 
into a solid (matrix) that can experience elastic deforma¬ 
tions in the form of local twists in the spirit of Einstein 
- de Haas effect. This leads to the time dependence of 
the orientation of the anisotropy axis. Rotations of the 
anisotropy axis are parametrized by the twist angle^l 

S(f{r,t) = X u{r,t), (2) 


where u{f, t) represents the displacement field of the crys¬ 
tal lattice surrounding the particle. Hamiltonian ([T|) is 
valid in the crystal frame. The corresponding Hamilto¬ 
nian in the laboratory frame can be obtained by substi¬ 
tuting n by the corresponding rotated anisotropy axis: 
the vector of local twists can be expressed as 5(j) = 0fc, 
with 6 — \54>\, and any rotated vector can be calculated 
through Rodrigues’ rotation formula: 


Vmt = vcosO + (k X v) sin0 -|- k{k ■ '(7)(1 — cos6*), 


V + 5(j) X V + -54> X {5(j> X v), 


( 3 ) 


where in the last step we have considered the regime of 
small deformations (|5^| 1) and an expansion up to 

second order in S(j). Therefore, the following expression 
for the Hamiltonian operator in the laboratory frame is 
obtained: 


Hlf — 'ft-cF — D{n ■ S,n X S} ■ 6(f — D[(^S(f x n) • ^ 

“ (4) 

= -DS^ + DiS^, Sy}6(l)^ - DlS^, S^}S(l)y 

- ^{Sz, Sa:}S(j)a;6(l)z - ^{S^,Sy}6(l)y6(j)z 

+ DSl {5^1 + 5cfl) - DSmi - 
+ D{Sx,Sy}S4>^S4>y, (5) 

where in the last step the anisotropy axis has been taken 
along the Z axis {n = e^) in the crystal frame. 

We consider two magnetic fields applied to the 
nanoparticle: the first one consists of a static field along 
the anisotropy axis in the crystal frame, whereas the sec¬ 
ond one is an AC field transversal to it. These fields are 
given by 


Ho = Hoe,, ( 6 ) 

Hi = K{t)ex, = Re , 

which add a Zeeman term to the total Hamiltonian, "Hz = 

-^hS- [Hq+Hi) . We will assume that both Hq and A are 
spatially homogeneous. From this point forward we also 
assume that the nanoparticle is located at the origin of 
coordinates, which translates into the local twists being 


evaluated at r = 0. Therefore, the Hamiltonian operator 
in the laboratory frame turns out to be 


"Hlf = -DSj - 'yfiS,Ho - jhSxhxit) 


( 7 ) 


- D 


/R3 


d3fj(3)(f^ {S,,Sx}5(j)y{f,t) 


- {Sz,Sy}54>x{f,t) A -{S„Sx}54>x{r,t)5(l),{r,t) 
+ ]^{Sz,Sy}5(l)y{r,t)5(l),{r,t) + 


- {Sx,Sy}S'px{r,t)S4>y{f,t) +Hph, 


where we have included the contribution H-ph to the total 
Hamiltonian from phonons. The classical spin is defined 

as S(t) = Tr[/3,^, where p is the density matrix oper¬ 
ator of the system, and the total energy of the system 
becomes: 


£: = Tr[p?iLF] (8) 

= fph - DZ,l{t) - 'yhZ,,{t)Ho - ^hY.x{t)hx{t) 





- 2'E,,{t)'Ey{t)d4>x{r,t) + T,,{t)T,x{t)6(j)x{r,t)d4>,{f,t) 

+ ^zit)T,y{t)S<py{r,t)S(l)z;{r,t) + Y.l{t)54>l{r,t) 

+ - E2(t)(5</.2(r,t) -h(5(/>2(r,t)) 


- 2Y.x{t)'Zy{t)5(t)x{r,t)5(j)y{r,t) 


Notice that in this treatment we are assuming that 
averaging over spin and phonon states can be done in¬ 
dependently, a legitimate assumption in the long-wave 
limit (|(5(/)| <C 1). As we are studying an AC-field-driven 
dynamical system, in the linear regime we have 'Zx^y A 
and u A. Therefore, we can linearize (up to first or¬ 
der terms in A) the dynamical equations arising from the 
total energy E due to the perturbative nature of this field. 

Landau-Lifshitz-Gilbert equation determines the dy¬ 
namical evolution of the classical spin: 


= (9) 


where //eff = —SE/SYl is the effective field stemming from 
the total energy of the system and rj is the dissipation 
constant. Eq. (|5]) and linearization leads to the following 
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system of equations: 

= 2DT.,T.y + -fhHoT.y + 2DT.l8(l)^{r = Q,t) 


rjh 


V _ V ^ 

^ dt dt 


( 10 ) 


_dYly 

dt 


= -2D^,i:,-jhHoJ:, + 2DEl68)yir = 6,t) 


+'yh^zhx{t) - 


rjh 


dT,j; dE^ 

" dt ^ dt 


’'"w = 


The third equation leads to the identity '£‘z{t) = Ez^o = 
cte up to second order in A. Furthermore, |E| ~ E^, p due 
to the X— and j/—components of the classical spin being 
perturbations. Introducing these dependences into the 
first two equations, we obtain the linear system of ODEs 

h-jA = hwFM^y + 2DElQS(l)oo{f=6,t)+r]h-^, 

= -fu^FM^x + 2DT,lQS(l)y{f=6A)-V^^^ 

+jhEz,fihx{t), (11) 

2DTjz 0 

where wfm = -^is the ferromagnetic reso- 

h 

nance frequency of the nanoparticle. 

The dynamical equation for the displacement field is 
given by 


d^Ua _ da^p 
'W ~ y dxp' 


( 12 ) 


where Oap = dEjSeap is the stress tensor, Cap = 
duajdxp is the strain tensor and p is the mass density. 
Notice that 8^ = (1/2) V x it, from which we straightfor¬ 
wardly deduce the identity 8(\)^ = {\/2)e-ypaeap. There¬ 
fore, the (linear) magnetic contribution to the stress ten¬ 
sor of the lattice is given by 

(iOcy/Sa (13) 

-b HE^ pd(^)(r) [8(j)x{rA)^xpa + 8(l)y(fA)>^vPa\ ■ 


III. BULK WAVES 


The dynamics of bulk acoustic waves are described by 
Eqs. (|12p . With account of expression (fT51) for the mag¬ 
netic contribution to the stress tensor, we obtain the fol¬ 


lowing dynamical equations for the bulk waves: 


9 Ux 2 v 72 _ ^ X' V ^ a(3)/'-a 

ot^ p OZ 


(14) 


H- \5^^\r)8(l)y{r,t) 


p dz - 

d Uy o __9 D „ d / —A 

^-CtV Uj, = --E,,oEj,—(r) 


DT.I. d r 


2^2 D 

-—-c,Vu, = -E,,o 


^ 8^^\f^8(j)xif,t) 

d 


p dz - 
d 






A 

dy L’ 


8^^\f)8(l)xif,t) 


A 

dx - 


8^^\f^8(l)yif,t) 


where we have set condition V • it = 0 for the displace¬ 
ment field. This stems from the fact that for nanoparti¬ 
cles, which are more rigid than their environment, only 
transverse phonons interact with the magnetic degrees 
of freedom, that is, local twists are transverse in nature. 
Performing the Eourier transform 


f{f,t) ^ fik,uj) = 


1 


(2^)= 


/r4 


dt d^f f{f,t)e 




on Eqs. (fTTl) we obtain the algebraic identities 
1 ikz 


(15) 


Ux{k,uj) = 


(27r)3/2p _ ^2 L 


— DUz,o'^xi^) 


(16) 


j{k,u;) = - 


Uz{k,uj) = 


ADY.l^F[8<l>y{r = 0,t)]{u:) 

ikz 


(27r)3/2pc2fc2 _^2 L 


DY.z,o81y(uj), 


+DY.l^F[8<l>x{?=Q,t)]{u^) 

D^z,o({ikx)'^x{^) + 


(27r)3/2p _ ^2 L 

{iky)ty{uj)) +DT,lQ({iky)F[8(j)x{f= 0,i)](a;) 


-(ikx)F[8(l)y{f = 0,t)](w)). 


With account of the identities (HU we can express the 

Fourier transforms of the twist angles in terms of E. 
Therefore, the above equations become 


Ux{k,uj) = 


j{k,uj) = 


1 


ikz 


(27r)3/2p p2^2 _ (^2 2 ’ 


1 


ikz 


z(k,uj) = - 


(27r)3/2p c2fc2 _ ^2 2 ’ 

1 h {ikx)a + iiky)b 


(27r)3/2p 2 Ap _ ^2 


( 17 ) 
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with 

a = —iuj'E. 


V 


2DS 


wfm 


z,0 


u ■ I I 2Z5S2 0 . 1 ^ 

0 = UjJIjx + 1 wfm-^-*W77 ) hy. 


Notice that the identities 
J^[6(j)^{6,t)]{uj) = 

T[S(l)y{0,t)](uj) = 


- iuJT] ) Sa; - 

(18) 


f i d^fc 

—> —# 

L 2(2^)3/2 

kyUz{k,uj) - kzUy{k,uj) 


i d^k 


(19) 

kzUxik,w) - kxUz(k,w 


/r3 2(2^)3/2 I 

hold, which combined with Eqs. (TlBll lead to 

h 


J^[54>x{f = 0,t)](a;) = 


(cc)6, (20) 


3(27r)2pc| 

T[5(l)v{r = d,t)\{uj) = - ^ S(a;)a. 

3(27r)2pc^^ 

The function S(a;) = (2/3)A:f£,+i7r(w/ct)^ stems from in¬ 
tegration over the 3D wave-vector space of the functions 

^2 I ^2 271- 

—5 I G {x,y}. The cut-off k^o = repre- 

C( fc2 — a;2 L 

sents the maximum modulus of a wave vector generated 
by the precession of an object of size L. These integrals 
can be calculated by means of contour integration in the 
complex plane, leading to the value {8tt/3c^)E{oj). 

In the linear regime the identities 

'Sx{oj) = Xxx{x})hx{uj), Sy(a;) = Xyx{xi)hx{uj), (21) 

hold. Taking into consideration the functional depen¬ 
dencies (1^01) of the magnetoelastic coupling terms, the 
dynamical equations for the spin susceptibility in the 
Fourier space become: 


ju;(l-|-C) Mi2 
—Mi2 ia;(l-|-C) 


\xx 

Xyx 


= -7^2,0 


0 

1 + C 


( 22 ) 


where 


We are interested in the regime |C| 1, within which 

functions a and b have the following functional depen¬ 
dence: 


(7S.,o) I — 


{u;FM-iujri)hxiuj)/K{uj), (26) 


/ 9 Tj'y \ 

b= (qSz.o) f— {-i^)hx{oj)lK{uj) 

Now we proceed to calculate the time and spatial de¬ 
pendences of the displacement field, that is, the profiles 
of the acoustic waves generated by the precessioiial mo¬ 
tion of a single spin. Introducing the above functional 
dependences for the spin susceptibilities into Eqs. dm, 
we obtain: 


c(d<) = 


1 


duje 


—iujt , 


Ak-r 


ikz 


(27r)7/2p 2 

ih f 

4(27r)3/2pc? 


a I d k e 

JR3 Cfk'^ — a;2 

dwe“*“*a ( —^ -ji[(a;/ct)r], 
CtJ r 


where in the last step we have used the planar wave ex¬ 
pansion 


e^k-f = dnY,^'ji[kr] ^ Yr{k)Yr{f) 


(27) 


z>o 


1——1 


for arbitrary position vector r, orthonormality of the 
spherical harmonics over the 3D sphere and contour inte¬ 
gration in the complex plane with respect to the variable 
k to calculate the integral over the 3D wave-vector space. 

The spherical Bessel function ji[x\ = ^^— xcosx 

termines the spatial dependence of the az-component of 
the displacement field. With account of the identity 
hx{oj) = (27r)3/2^^(i) — for a circular AC field, we fi¬ 

nally obtain after integration over the frequency domain: 


Ux(.r,t) = j 


Sirpc] 


Ct 


(28) 


with A G R. We can proceed analogously with the y and 
z-components of the displacement field: 


, , 2DY,,o , , 2DE2 S(a;) %(d0 = 

Mi2 = (a;FM-*a;?7)(l-fC)-^C, C w = ! ^ ■ 

h 6[27T)^pC^ 

(23) Uz{f,t)=- 

This system of equations can be solved using standard 
methods, so that we obtain the following expressions for 
the spin susceptibility 




Sirpc'^ V 


) ji[(C/ci)r], 


'yhEzfiA f 2DT,z^o 
St: pc. 


b 




X I ^ ) ji[{Act)r\ 


\xx 

Xyx 


where 


7^z,o(l + C) f —Mi2 

Kioj) \ *a;(l -|- C,) 


K{u)=A{l + 0-Ml2. 


(24) 


(25) 


with 


A[^, i\ = Re 
i?[^, t] = Re 


OZFM - 

tK{0 


(29) 


( 30 ) 


-e - 


KiO 
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where K{uj) = — (wfm — +a;^77^) + 

2iu:ujYM'n- resonance we have K{^ wfm) — 2 ia;pj^? 7 , 
where condition ry <C 1 has been taken into consideration. 
Therefore, the expressions for the functions A and B close 
to resonance are: 


A[^ ~ iOFM,t] = 7; -[77sin(a;FMi) - cos(a;FM0]> 

ZujfmV 

B[^ — wfm, t] = - -sin(wFM0- (31) 

ZUJfmV 


IV. SURFACE WAVES 


The geometry of the problem is the following: the sur¬ 
face of the solid is located at the XY plane and its bulk 
extends over the Z < 0 space. The mechanical stress 
tensor is given by: 

E 

^XX ~ ^2^ _|_ — 2cr) ~ ^)'^XX + ^{Uyy + Uzz)\ , (32) 

m ^ 

O'yy = — 2cr) ~ ’^{'^XX + Uzz)] , 

^ZZ ~ ^ 2 ^ _j_ — 2cr) ~ ^)'^ZZ + 0'(uxx + Uyy)\ , 

_m E E _ E 

+ (l + cr)''""’ + 

where = \{dpUa + d^up), E is the Young modulus 
and cr is the Poisson ratio. Notice that the transverse and 
longitudinal speeds of sound are given by cf = E/2p{l -|- 
cr) and cf = E{1 — a)/p{l + cr)(l — 2a), respectively. 

Boundary conditions at the surface are given by 
o^rij = 0, where n is the normal vector to the surface. 
In our case rij = 6j^z, so that tr™ = cr™ = cr™ = 0 at 
z = 0. The total displacement field is a superposition of 
transversal and longitudinal modes, with the transversal 
mode given by the bulk solutions (E51) and The ex¬ 
pected functional form for the longitudinal surface wave 
is given by 


(27r)V^ Jr3 

(33) 

where R = {x,y), q = {qx,qy) is the 2D wave vector 
and OL{q) is the inverse of the penetration length. The 
longitudinal mode satisfies the condition V x ui = 0, 
which in the Fourier space turns into the set of equations 

iqyiiz - avy = 0, (34) 

iqxVz - avx = 0, 

qxVy - qyVx = 0 . 


boundary conditions in the Fourier space become 

1 


(27r)i/2 

f a dk 


/ dkz [{ikz)ux,t + {iqx)uz,t] + 2aVx = 0, (35) 
Jr 


dkz [{ikz)uyp -I- {iqy)uz,t] + “^avy = 0, 

\jdkz)'Uzp “t” (1 2A ) {{'iqx)‘^x,t d~ {iqy)uypy^ 


Jr (2^)i/2 

+ \c^ + (2A^ — \)(f'\vz = 0, 

where A = ct/ci is a parameter of the solid. Computation 
of the integrals leads to the following results: 


( 2 : 


.)l/2 — ^27r)l/2 

—h iqxtt + iqyb 


1 


(2^)i/2 7 r 

1 f 


dkz Uz,t = 
dkz {ikz 
dkz {ikz 


(27r)i/2 


'R 


Sirpc^ at{q) 

ha 

■ 


at - 

hb 


Sttpc? 

at - 

z,t = 0, 


ct)^] and kiD 


2fciD 


2k 


ID 


cut-off wave vector along the Z axis. With account of 
these expressions, the third equation of the set ((35ll be¬ 
comes \a^ + (2A^ — l)q'^]vz = 0, which leads to the follow¬ 
ing dependence for the inverse of the penetration length: 


a{q)= [l-2Xp^%. 


(37) 


Combining the first two equations of (13511 with account 
of Eqs. (1341) yields the identity 


Vziq,uj) = 


-q 


2q2 


dkz 


Uz,t 


iqy 


JR 

iqx 


'R 


R 


IbTrpCj 


2kiD , q 

at -1- 

TT at 


dkz{ikz)ux,t 

^'' iqxO- + iqyb 


dkz{ikz)uy,t 


(38) 


where in the last step we have used Eqs. (IMll . The func¬ 
tional dependence of Vz on the surface can be obtained 
by inverse Fourier transformation: 


,{Ra = 0,t) = 


f 


(27r)3/2 


'R3 


da;dVwz(9,a;)e*('?-^-“*) 


( 39 ) 


ih r 

8(27r)5/2pc? 


d^qm{q,uj)^E^- 


+ 


ih 


8(27r)3/2pc' 


dw I 


t JR 


/r2 9 

d^q m{q,uj)Jje 


/R2 


q 


With account of all these functional dependences, the 
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where m{q,uj) = at{q) - 

77 

the Jacobi-Anger expansion 


—;rw. With account of 

at[q) 


^ ^ ent^MqR) cos [n{9 - 6 r)] , 

n=0 

(40) 

where Cq = 1, e„ = 2, Vn > 1, J„(a;) is the Bessel func¬ 
tion of order n and 9r is the angle between R and the X 
axis, and with account of the orthogonality of the Fourier 
basis {1} U {sin(n0), cos(n0)}„>i, we have that 



47riG[i?, w] 


cos 9 r 1 
sin 0r J ’ 


(41) 


where G[R,uj] = 


dg m{q,uj)Ji{qR) and q 2 D = 


^{2Tr) jL is the corresponding two-dimensional cut-off 
wave vector. Therefore, Eq. becomes 


Vz{R,z = 0,t) = 


SirpCf V ^ 




(42) 


Generation of these acoustic waves by a precessional 
nanoparticle can be detected experimentally by means 
of the TEM set-up. The current spatial resolution is 
about 0.01 A. Potential candidates to host an observ¬ 
able acoustic deformation of the substrate should com¬ 
bine high (uniaxial) anisotropic nanoparticles and a soft 
substrate. In the following, we focus on the case of an 
FePt nanoparticle deposited on a gold substrate. 

As an elastic medium, gold is characterized by the val¬ 
ues p = 19.3 g/cvD? for the mass density and ct = 120000 
cm/s for the transversal speed of sound. As a conse¬ 
quence of the crystalline ordering of the face-centered 
tetragonal structure, FePt nanoparticles in the LIq phase 
show a very high uniaxial magnetocrystalline anisotropy. 
The corresponding lattice constants are a = 3.85 A and 
c = 3.71 A. Within a tetragonal cell there are 4 Pt atoms 
(at the center of the lateral faces) and 10 Fe atoms (2 
at the center of the top and bottom faces and 8 at the 
vertexes of the cell). With account of the spin S' = 5/2 
for these atoms, the spin density of the Llo-phase FePt 
nanoparticles is = 1-82 • 10^^ cm“^. The ferro¬ 

magnetic resonance frequency of these nanoparticles is 
wfm = 5.63-10^° Hz. Therefore, assuming a cubic geom¬ 
etry of the spin cluster we have where 

L is the size of the nanoparticle. 


f- 

\R 


M.t] + ^B[U]) ImG[S,e] 


where 


G[^,t] =Re 
D[^, t] = Re 


i^rj - CUFM -i^t 
" . 


(43) 


- 




KiO 

At resonance, these functions become 

G[^ ~ WFM,^] = 7^-[sin(a;FMi) + P cos(a;FMi)] 

ZwfmV 

D[^ — wfm, t] = - -cos(a;FMi)- 

2wfmP 


(44) 



V. DISCUSSION 

We have computed acoustic waves generated by the 
spin precession. Our results are parameter-free in a 
sense that they do not contain any unknown parame¬ 
ters of spin-phonon interaction. The computed ampli¬ 
tudes of the acoustic waves have been expressed in terms 
of parameters that can be independently obtained from 
macroscopic experiments, such as FMR frequency, speed 
of sound, etc. The effect comes from the elastic twists 
generated by torque associated with the spin precession. 
It is a microscopic counterpart of the Einstein - de Haas 
effect. Its observation would provide an interesting test 
of theoretical concepts. 


Figure 1: (Color online) Functional dependence of the maxi¬ 
mum amplitude of both bulk (black) and surface (red) acous¬ 
tic waves on the size L of the nanoparticle when measured 
at the distance R — 2L. The blue dotted line represents the 
resolution 0.01 A of the STM technique. 

As a criterion to determine whether the effect is mea¬ 
surable we estimate the maximum amplitude of the in¬ 
duced acoustic wave at the distance from the particle 
that equals its radius and compare it to spatial resolu¬ 
tion of the STM technique. This basically suggests that 
there will be a measurable effect in the immediate vicinity 
of the particle within the nanoscale range. Fig. 1 shows 
the dependences of the maximum amplitude of both bulk 
and surface contributions to the substrate elastic defer- 
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mation on the size of the nanoparticle when measured at 
the distance R = 2L, 


KmLi = ' |ji[(a.FM/Q)2L]|, (45) 

167r77pcf \ ct J 

l<”ax“l = ^^^V[ReG[2L,a;FM]]^ + [Im G[2L, a.FM]]^ 


where we have approximated 


2i4E2^o 

h 


wfm {Ho ^ 100 


Oe) and considered rj = 0.001. Notice that for sizes in the 
ballpark of L = 135 nm both bulk and surface acoustic 
waves show a maximum amplitude bigger than the spa¬ 
tial resolution of the STM technique, so that these elastic 
deformations could be experimentally detected. The de¬ 
tection of only induced surface acoustic wave could be 
performed by either considering bigger sizes (in the ball¬ 
park of L = 180 nm, where the maximum amplitude of 
the bulk wave is strongly suppressed) or using nanoparti¬ 
cles with smaller ferromagnetic resonance frequency. The 


latter is due to the explicit quadratic dependence of Eq. 
(1451) on this frequency whereas the maximum amplitude 
of the surface contribution is less sensitive to it. 


Expression (E51) can be approximated by C(‘^fm) — 

^FM z o (^fm) FePt nanoparticles we 

3(27r)2pc^ 

have Id < 7 • 10“'^ within the range of sizes consid¬ 
ered above. Therefore, the results obtained are consistent 
with the assumption |d “C 1 made in Sec. 3. Note that 
particles detectable by means of TEM via acoustic waves 
generated by the spin precession would be much smaller 
than samples required for the detection of the ESR signal 
in conventional macroscopic measurements^^. 
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